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The Advanced LIGO and Advanced Virgo ground-based gravitational-wave (GW) detectors are 
projected to come online 2015-2016, reaching a final sensitivity sufficient to observe dozens of binary 
neutron star mergers per year by 2018. We present a fully-automated, targeted search strategy for 
prompt gamma-ray counterparts in offline Fermi-GBM data. The multi-detector method makes use 
of a detailed model response of the instrument, and benefits from time and sky location information 
derived from the gravitational-wave signal. 


I. INTRODUCTION 

Compact binary coalescence (CBC), such as the 
merger of two neutron stars (NS) or black holes (BH), 
remains the most highly anticipated gravitational- 
wave signal for ground-based gravitational-wave de- 
tectors. The second-generation Advanced LIGO 
(4km-baseline interferometers in Hanford, WA and 
Livingston, LA) [1] and Advanced Virgo (3km inter- 
ferometer in Cascina, Italy) [2] detectors are projected 
to come online 2015-2016, and reach a final sensi- 
tivity sufficient to observe dozens of NS-NS mergers 
per year by 2018. Together with the 600m GEO-HF 
detector [3], they will form a word- wide network of 
gravitational-wave interferometric detectors. The net- 
work will be joined by the Japanese 3km cryogenic 
KAGRA detector [4] and a proposed third LIGO-India 
detector around 2018-2020, increasing overall sensi- 
tivity and sky-coverage, while also improving sky lo- 
calization and waveform reconstruction [5]. 

Gravitational-waves have yet to be directly ob- 
served. Our confidence in the existence of com- 
pact binary coalescence events comes primarily from 
the discovery of a small number of galactic pulsars 
which appear to be in close binary systems with an- 
other neutron star. The first and most famous of 
these systems contains the Hulse-Taylor pulsar PSR 
B1316+16. It has demonstrated over many decades an 
orbital decay consistent with loss of energy to gravi- 
tational waves [6]. The number and inferred lifetime 
of these systems can be used to obtain an estimate of 
about one NS-NS merger event per Mpc 3 per million 
years [7], with up to two orders of magnitude uncer- 
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tainty in rate due largely to our limited knowledge 
of the pulsar luminosity function and limited statis- 
tics. The merger rate translates into an estimate of 
~0.02 detectable NS-NS merger events per year for 
the initial LIGO-Virgo detector network, which op- 
erated between 2005-2010, and ~40 per year for the 
advanced detectors once they reach design sensitivity. 
Although their gravitational radiation is stronger, the 
merger rates of NS-BH binaries is more uncertain as 
we have not observed any NS-BH binary systems, and 
have generally poor knowledge of the black hole mass 
distribution. 

Gamma-ray bursts (GRB) are flashes of gamma 
rays observed approximately once per day. Their 
isotropic distribution in the sky was the first evi- 
dence of an extra-galactic origin, and indicated that 
they were extremely energetic events. The duration 
of prompt gamma-ray emission shows a bi-modal dis- 
tribution which naturally groups GRBs into two cat- 
egories [8]. Most long GRBs emit their prompt ra- 
diation over timescales >2 seconds, and as much as 
hundreds of seconds. They have been associated with 
young stellar populations and the collapse of rapidly 
rotating massive stars [9]. Short GRBs (sGRB), with 
prompt emission typically less than Is and a generally 
harder spectrum, are found in both old and new stellar 
populations. Mergers of two neutron stars, or of neu- 
tron star/black hole systems, are thought thought to 
be a major contribution to the sGRB population [10]. 
It is this favored progenitor model which makes short 
GRBs and associated afterglow emission a promising 
counterpart to gravitational-wave observations. 

Since 2005, the Swift satellite has revolutionized 
our understanding of short GRBs by the rapid ob- 
servation of x-ray afterglows, providing the first lo- 
calization, host identification, and red-shift informa- 
tion [11], The beaming angle for short GRBs is highly 
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uncertain, although limited observations of jet breaks 
in some afterglows imply half-opening angles of Oj ~ 
3-14 degrees [12, 13]. The absence of an observable 
jet break sets a lower limit on the opening angle which 
is generally weak (due to limits in sensitivity), though 
in the case of GRB 050724A, late-time Chandra ob- 
servations were able to constrain Oj > 25° [14], 

The observed spatial density of sGRB’s and lim- 
its on beaming angle result in a NS-NS merger event 
rate roughly consistent with that derived from galactic 
binary pulsar measurements. Although the beaming 
factor of ~ Oj / 2 means we believe most merger events 
seen by the advanced GW detectors will not be accom- 
panied by a standard gamma-ray burst, this is some- 
what compensated by the fact that the ones that are 
beamed toward us have stronger gravitational-wave 
emission. Current estimates for coincident GW-sGRB 
observation for advanced LIGO-Virgo are a few per 
year assuming a NS-NS progenitor model [15, 16]. The 
rate increases by a factor of 8 if all observed short 
GRBs are instead due to NS-BH (10 M 0 ) mergers 
which are detectable in gravitational-waves to twice 
the distance. 

In addition to the jet-driven burst and afterglow, 
other EM emission associated with a compact merger 
can be a promising channel for GW-EM coincidence, 
particularly if the EM radiation is less-beamed or even 
isotropic. A few short GRBs (~10%) have shown 
clear evidence of high-energy flares which precede the 
primary burst by 1 TO seconds, and possibly up to 
100s [17]. The precursors can be interpreted as evi- 
dence of some activity during or before merger, such 
as the resonant shattering of NS crusts [18], which 
could radiate isotropically. Thus it will be interesting 
to search for weak non-standard EM emission accom- 
panying all nearby NS-NS mergers seen in GWs, while 
we expect only a small fraction to be oriented in out 
line-of-sight for a standard jet-driven sGRB. 

The Gamma-ray Burst Monitor (GBM) [19] aboard 
the Fermi spacecraft measures photon rates from 8 
keV-40 MeV. The instrument consists of 12 semi- 
directional Nal scintillation detectors and 2 BGO scin- 
tillation detectors which cover the entire sky not oc- 
cluded by the Earth (about 65%). The lower-energy 
Nal detectors have an approximately cos 0 response 
relative to angle of incidence, and relative rates across 
detectors are used to reconstruct the source location 
to a few degrees. The BGO detectors are much less 
directional, and are used to detect and resolve the 
higher energy spectrum above ~200 MeV. 

GBM produces on-board triggers for gamma-ray 
burst events by looking for multi-detector rate ex- 
cess over background across various energy bands and 
timescales. In the case of a trigger, individual pho- 
ton information is sent to the ground and the event 
is publicly reported. Those events which have been 
confirmed as GRBs have already been studied in co- 
incidence with LIGO-Virgo data [20-22]. So far, no 


gravitational-wave counterparts to triggered GRBs 
have been identified, which was not unexpected given 
the limited reach of the first generation instruments. 

In addition to the triggered events, survey data is 
available which records binned photon counts over all 
time. In this proposed offline analysis of short tran- 
sients, we consider the CTIME daily data, which con- 
tains counts binned at 0.256s over 8 energy channels 
for each detector. A new GBM data product (contin- 
uous TTE) was implemented in late 2012. It provides 
continuous data on individual photons with 2/xs and 
128 energy channel resolution, which will further en- 
hance offline sensitivity to particularly short bursts. 

This paper discusses the possibility of using these 
offline GBM data products to follow-up gravitational- 
wave candidates in the advanced LIGO-Virgo era. In 
section II we describe the characteristics of a trigger 
provided by gravitational-wave data. Sections III and 
IV develop a likelihood-ratio based procedure for ana- 
lyzing available GBM offline data about the time and 
sky location provided by the gravitational waves. Fi- 
nally, in section V we demonstrate the performance 
of the search algorithm on background times and a 
sample of known sGRB’s. 


II. GRAVITATIONAL- WAVE TRIGGER 

Searches for NS-NS and NS-BH coalescence in 
gravitational-wave data typically use matched filter- 
ing of the model waveform [23, 24]. The gravita- 
tional waves from CBC are characterized by a chirp 
of monotonically increasing frequency and amplitude 
as the binary inspirals under approximately adiabatic 
orbital decay. During this progression, kilometer-scale 
interferometers are sensitive to the inspiral phase of a 
stellar-mass coalescence until just prior to merger, but 
before any tidal effects on a NS become important. 

The waveform in this regime is well-modeled us- 
ing post-Newtonian techniques, and can be consid- 
ered a standard candle encoding orbital parameters 
of the system. Mass, spin, and coalescence time are 
accessible with the GW projection onto a single de- 
tector, while distance, inclination, and sky location 
degeneracies can be disentangled to varying degrees 
of success using coherent data from multiple detec- 
tors [25]. A typical gravitational- wave detection in 
the early advanced detector era may have moderate 
signal-to-noise ~8 in two or more detectors, estimated 
merger time to within milliseconds, and sky localiza- 
tion to 'MOO square degrees depending on detector 
network configuration [26]. Figure 1 shows a recon- 
structed 3-detector localization for a simulated signal 
near detection threshold. 
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simulated NS/NS merger localization 



FIG. 1: Example sky localization of simulated NS-NS 
merger seen in modeled gravitational-wave data with noise 
similar in spectrum to that expected from the advanced 
LIGO-Virgo GW network. The signal-to-noise of the GW 
signal in each detector is 8.7 (HI) 7.2 (LI) and 3.1 (VI), 
and the sky location is determined through Bayesian in- 
ference using a nested sampling technique [25] applied to 
data from all detectors. The elongated shape beyond In- 
is due to a small degenerate timing uncertainty in VI. 


III. COHERENT ANALYSIS OF GBM DATA 

In this section, we develop a procedure to coher- 
ently search GBM detector data for modeled events. 
The basic idea is that by processing multiple detec- 
tor data coherently, we can obtain a greater sensi- 
tivity than when considering one detector at a time. 
Greater computational resources available offline (vs. 
on-board) also allow for more careful background es- 
timation to be done. For this analysis, we can relax 
to some extent the strict 2-detector coincidence re- 
quirement used to veto spurious events on-board as 
the gravitational-wave trigger means much less time 
and sky area is considered. 


where T ~ls, we estimate the background using a 
polynomial fit to local data from [— 10T, +10T] (min- 
imum ±5s), excluding time [— 3T/2, +5T/2] around 
the foreground interval to avoid bias from an on-source 
excess. An example of the foreground and background 
intervals about a strong prompt signal is shown in fig- 
ure 2. The polynomial degree is determined by the 
interval length to account for more complicated back- 
ground variability over longer intervals. It ranges from 
2 (minimum) to 1+0.5 log 2 T. The quality of the fit 
is determined using a y 2 statistic applied to the data, 
re-binned at T/4 seconds for T > Is. If the y 2 per 
degree-of- freedom is over 2, or if any of the N indi- 
vidual data points has y 2 > 4 In A, it is assumed that 
the polynomial could not adequately model the local 
background variability, and the fit is redone over the 
smaller interval [— 5T, +5T]. If the fit continues to fail 
the y 2 test at the looser requirements of 3 and 6 In N 
respectively, the background estimate is marked un- 
reliable and not used in further computations for that 
particular foreground interval. 

High-energy cosmic rays striking a Nal crystal can 
result in long-lived phosphorescent light emission. 
The detector may interpret this is a rapid series of 
events, creating a short-lived jump in rates for one or 
multiple channels, and severely distorting the back- 
ground fit if not accounted for. They are identified 
with a simple procedure that compares the counts in 
each 0.256s bin against the mean rate estimated from 
four neighboring bins. If an excess is detected to 5er or 
more, that measurement and the immediately neigh- 
boring bin on either side is removed from the back- 
ground fit for all channels of the affected detector. 
A much stricter requirement is used when rejecting 
cosmic-ray events in the foreground interval. In that 
case, a 0.256s measurement must have signal-to-noise 
over 500 /\/T for it to be excised, which is much larger 
than expected for true GRBs. 


B. Likelihood-ratio statistic 


A. GBM background estimation 

Each detector is subject to a substantial time- 
varying background from bright high-energy sources 
that come in and out of the wide field of view, as 
well as location-dependent particle and Earth atmo- 
spheric effects. This background must be estimated 
and subtracted out to look for any prompt excess. 
Methods in use include the local averaging of pre- 
vious data done on-board, smooth spline-fits with a 
high-frequency cutoff, direct tracking and modeling of 
the dominant sources, and averaging rates from pre- 
vious orbits with similar orientations [27 29]. In this 
analysis where we are interested in the background es- 
timate for a short foreground interval [— T/2, +T/2] 


A likelihood ratio combines information about 
sources and noise into a single variable. It is defined 
as the probability of measuring the observed data, d, 
in the presence of a particular true signal Hi (source 
amplitude s > 0) divided by the probability of mea- 
suring the observed data in noise alone Hq (s = 0). 


A(d) 


P{d\Hi) 

P{d\Ho) 


( 1 ) 


When signal parameters such as light-curve, spec- 
trum, amplitude s and sky-location a, 6 are unknown, 
one can either marginalize over the unknown param- 
eters, or take the maximum likelihood over the range 
to obtain best-fit values. 
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FIG. 2: Detector and channel-dependent polynomial back- 
ground fitting applied to GRB 090510A. Four of the 8 x 14 
channel and detector combinations are shown. The green 
represents the highest likelihood foreground interval, and 
the background contribution is estimated using a polyno- 
mial fit to the region (cropped) shown in blue. Red dots 
mark the average flux in the foreground interval. 


For i measurements of sufficiently binned, uncorre- 
lated Gaussian data, 


P(d i \H 1 ) = l 


1 

V2n(Td, 

1 


P(d t \H 0 ) = l[^= 

v Z7T< 


exp 


exp 



( 2 ) 

( 3 ) 


where we have used di = di — ( rii ) to represent the 
background-subtracted measurements (e.g. red dots 
minus blue curve in figure 2), a ni and a ^ for the stan- 
dard deviation of the background and expected data 
(background+signal) , ri for the location/spectrum- 
dependent instrumental response, and s, a single in- 
trinsic source amplitude scaling factor at the Earth. 
Maximizing the likelihood ratio is the same as maxi- 
mizing the log-likelihood ratio C = In A, 




In - 


di (di - ns ) 2 


a di 


2 cr 2 


2 er? 


( 4 ) 


When the first two terms are fixed, maximizing the 
log-likelihood ratio is equivalent to minimizing the the 
third term which we recognize as a x 2 fit to the two 
free sky-location parameters in ri and the single am- 
plitude parameter s. The shape of the likelihood func- 
tion over source amplitude s is the product of (almost) 
Gaussians centered about the best-estimate of s in 
each measurement, and the total likelihood is scaled 
by the sum of the squared signal-to-noise of all mea- 
surements. 


The dependence of response factors n on sky loca- 
tion is complicated, so the likelihood ratio is calcu- 
lated over a sample grid of all possible locations. As- 
suming a single location, the remaining free param- 
eter is the source amplitude s. The variance in the 
background-subtracted detector data includes both 
background and source contributions, 

Vdi = + OS + cr 2 .s 2 (s > 0) (5) 

with it 2 , representing Gaussian-modeled systematic 
uncertainty in the instrumental response. Source 
terms are only included for physical s > 0 else their 
contribution is zero. The background contributes 
Poisson error, as well as any systematic variance er 2 
from poor background fitting which is also assumed 
to be Gaussian, 

vli = ( n i) + °i- ( 6 ) 

We find Sbest which maximizes C by setting the deriva- 
tive dC/ds to zero. If = cr ni , which happens when 
rs -C cr 2 , v / 2^max reduces to a coherent SNR (sum 
data using weights g/ct 2 .). If can be assumed 
constant, the solution for Sbest is found analytically, 


E; r idi/°di 

Ei r?/a 2 di 


( 7 ) 


which is just an appropriately inverse-noise weighted 
sum of the individual estimates of s from each mea- 
surement. Although depends on s, we can make 
the practical approximation cr^. rs ma x((rii},di) + er 2 
as an initial guess. 

To find the true maximum, we solve for dL/ds = 0 
using Newton’s method, beginning with the value of 
the first derivative at our initial guess soi an d using the 
analytic second derivative to refine the measurement: 
Sbest ~ so — (dC / ds) / (d 2 C / ds 2 ) . This calculation is 
both fast and easily vectorized. One initial guess using 
equation 7 followed by a couple iterations of Newton’s 
method generally provides an excellent approximation 
for Sbest ■ 

So far we have been calculating the probability of 
a signal assuming a specific source amplitude s. To 
consider all possible source amplitudes we need to in- 
tegrate the likelihood P(d\s) (equation 2) over a prior 
on s, 


P(d) 


P(d\s)P(s)ds 


( 8 ) 


For a given set of detector data d , the likelihood P(d\s) 
over s is almost the product of individual Gaussian 
distributions (not quite Gaussian because ad depends 
on s). The product of Gaussian distributions with 
mean values fii and standard deviations <Ti is itself 
Gaussian with mean and variance, 


_ E in/vi 

Mprod — ^ l jfj 2 ’ 


a 


2 

prod 


l 

EW 


(9) 
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In this case, pi = d{ /r,; and cr, = a di/^i- The mean 
value //prod is the same as our initial guess for s at 
maximum likelihood (equation 7), but we can assume 
to have a more accurate maximum-likelihood location 
Sbest from the numerical procedure outlined above. 
The estimate for the variance of £ over s is, 


Erf !a\ 


af. evaluated at s b est (10) 


For a flat prior P(s) = 1, we can integrate equa- 
tion 8 by simply considering the area of a Gaussian 
with peak value P(d|sb es t) and variance erj. Other 
choices for a prior can be represented by a power-law 
distribution P(s) oc s A spatially homogeneous 
population, suitable for nearby sources, would assume 
/3 = 5/2, while the empirical distribution for observed 
GRB amplitudes follows a power-law decay closer to 
/3 « 1.8, reflecting cosmological effects. If instead we 
are looking for signals from a particular host galaxy 
at known location and distance, we want to use an 
intrinsic source luminosity distribution. A convenient 
option is to use a scale-free prior with fixed /? = 1, 
so that our choice of form for an amplitude prior at 
the Earth does not translate back into a luminosity 
distribution that varies with distance. 

One difficulty with any power-law prior is that it 
diverges for s — > 0. In a reasonable scenario of SNR 
>a few, the integral over s will consist of two distinct 
contributions. The first is a Gaussian component with 
p = Sbest and variance crj, scaled by the prior P(sbest)- 
The second component is an infinite contribution from 
the divergent prior at s = 0 with little contribution 
afterward due to suppression from the Gaussian tail. 
The infinite contribution represents the certainty of a 
signal of arbitrarily small amplitude to be present in 
the data, regardless of the ability of the data to resolve 
it. For moderate SNR, it’s very easy to isolate only 
the Gaussian contribution by placing a small cut on 
amplitude, truncating the likelihood for s < s cu t- Un- 
der the approximation that P(s) varies slowly over the 
width of the Gaussian, the log-likelihood marginalized 
over all s becomes, 


£(d) = -/3 In Sbest +ln dc + £(d\s hest ) (11) 


up to common additive constants that do not depend 
on the data d. 

The approximation for the marginalized likelihood 
becomes problematic for small Sbest where the as- 
sumption that P(s) is constant over the range of the 
Gaussian breaks down, and the Gaussian distribution 
can no longer be isolated form the divergence at small 
s. We enforce a finite and well-behaved prior by mul- 
tiplying by a prefactor, 


P(s) 


X _ g — (s/T^-cU 


( 12 ) 


so that P(s) reaches a maximum constant value of 
(l/'yac)^ for small s. The tunable parameter 7 sets 
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FIG. 3: Signal-to-noise expected from GRBs with a nor- 
mal spectrum, and coming from various positions across 
the GBM field-of-view. The hypothetical signal lasts 
0.512s and is normalized to 1.0 photons/s/cm 2 in the 50- 
300 keV band, while background rates and Earth position 
are selected from 10 seconds prior to GRB 090305A. The 
instrumental response includes contributions from atmo- 
spheric scattering in the 50-300 keV band. The first map 
shows the signal-to-noise from the single best detector over 
50-300 keV. The second plot shows the SNR expected from 
a coherent analysis of the full CTIME data. The coher- 
ent statistic must be calculated for all directions, which 
means it is subject to a trials factor equal to the number 
of independent sky locations. 


the number of standard deviations at which the prior 
begins to plateau, and we use 7 = 2.5. This allows 
us to use the approximation that P(s) is reasonably 
constant over a range of <j£ for any s > 0. The only 
remaining correction is to account for clipping of the 
Gaussian for non-physical s < 0, which can be rep- 
resented by the error function. The final approxima- 
tion for the amplitude-marginalized log-likelihood be- 
comes, 


£{d) = lno£ + In 


1 +Erf 


'Sbest 


1 In 

X _ g-lsbest/TO’c)' 3 

l-ifi 

In (7 <t £ ) 


j 1 ' -£(d|.Sbest) 

/31n Sbest Sbest A 0 
Sbest A 0 


(13) 


which contains factors form the Gaussian width, frac- 
tional overlap with s > 0, maximum likelihood at 
Sbest) and scaling from P(s) respectively. Finally we 
are free to calibrate the log-likelihood by subtracting 
the expected £{d) calculated for no signal at a ref- 
erence sensitivity: £ re f = — /31n 7+ (! - /3) lucre, ref • 
0£ represents the source amplitude required for a 
ler excess in the combined data, and is around 0.05 
photons/s/cm 2 x (T / Is) -1 / 2 [50-300 keV] for typi- 
cal source spectra and reference background levels. 
Figure 3 shows the coherent signal-to-noise expected 
from all detectors for a 0.512s-long event with nor- 
mal GRB spectrum and constant amplitude of 1.0 
photons/s/cm 2 , and compares it to the SNR expected 
from the most favorably-oriented individual detector 
alone in the 50-300 keV band. 
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IV. PERFORMING THE FOLLOW-UP 

For a single gravitational-wave trigger, we search 
a standard sGRB accretion timescale of [0, 5s] rel- 
ative to the time of coalescence for prompt flux ex- 
cess in GBM between 0.256 and 2s long (the lower 
limit of 0.256s set by the CTIME archival resolution 
does not apply to continuous TTE data). We also 
search a prior interval [-30s, Os] for possible precur- 
sor bursts between 0.256s and 2s. Finally we may 
include bursts between 2s to 32s in an extended inter- 
val [ -30s, 300s] to search for possible longer-duration 
emission. While emission outside of the standard ac- 
cretion timescale can be considered speculative, it’s 
worth looking for given that any events detectable in 
gravitational- waves will be closer than sGRB’s with 
known red-shift to date, making them good candidates 
to search for weak exotic and possibly less-beamed 
emission. To appropriately tile the search in time and 
duration T, we use rectangular windows with T spaced 
by powers of two (0.256s, 0.512s, Is, etc.). Their cen- 
tral times are sampled along the search interval in 
units of T /4 to provide an even mismatch in signal- 
to-noise across search windows. 

The likelihood-ratio statistic described in the pre- 
vious section is calculated for each foreground inter- 
val based on the background-subtracted flux measure- 
ments di in each of the 8 x 14 channel and detec- 
tor combinations. In addition, the instrumental re- 
sponse Ti (and thus the likelihood-ratio) also depend 
on source sky location and spectrum. To efficiently 
calculate the response over a large area of sky, we use 
precomputed all-sky response look-up tables originally 
generated for offline localization [30]. The likelihood 
ratio as a function of sky position provides a proba- 
bility distribution over the sky of a GBM signal (as in 
figure 5). This can be coincided with the GW-derived 
skymap (figure 1) by direct multiplication, in effect 
using the GW skymap as a prior. 

Finally, we marginalize over sky location and rep- 
resentative source spectra to get a ranking of events 
characterized only by their foreground time intervals. 
A unique list of non-overlapping events is constructed 
by beginning with the highest-ranked event and re- 
moving from the list any lower-ranked events with 
foreground windows that overlap it; then taking the 
next surviving highest ranked event, etc., until the list 
is exhausted. 


V. TEST ON SWIFT SHORT GRBS 

We can test the offline analysis on short GRBs trig- 
gered by Swift and observable by Fermi-GBM. The 
Swift GRBs are particularly useful as their compara- 
tively good localization can be used to test systematic 
errors introduced through the GBM model response. 



FIG. 4: The maximum likelihood coherent flux excess 
found by searching time around Swift short-GRBs within 
the GBM field of view using foreground intervals between 
0.256 and 2s. GRBs are placed at 0.01 Hz in anticipation of 
a GW-GBM coincidence search window of <100s. An es- 
timated cumulative background rate distribution is found 
by taking the non-overlapping all-sky events from a con- 
tinuous scan of approximately one-day in total of nearby 
data. The dotted background distribution shows the ef- 
fect of a simple cut to remove particle events (the cut does 
not remove any GRBs). GBM excesses corresponding to 
GRB 090815C and 110112A are buried within the back- 
ground distribution, and are not clearly associable with 
their GRBs. The ability of the method to identify the 
next two weak, un-triggered events GRBs 090305A (see 
also figure 5) and 120403A depends on ongoing strategies 
to further reject non-Gaussian outliers. The remaining 
Swift sGRB’s in the sample have a corresponding GBM 
on-board trigger. 


We calculate the likelihood ratio for all sky positions 
for GBM data local to Swift sGRB measurements (fig- 
ure 4). The all-sky likelihood ratio is proportional to 
the probability that the source originates from a par- 
ticular location on the sky, and can be used to produce 
a GBM skymap. 

About 50% of Swift sGRBs are within GBM’s field- 
of-view (the 65% of the sky not occulted by the Earth) 
and occur during an interval of time when Fermi is 
outside of the South-Atlantic Anomaly (SAA) and 
operational. Most of these observable Swift sGRBs 
also triggered GBM on-board. GRB 081024A trig- 
gered on-board, but was too close to the interruption 
of data taking during passage through the SAA for our 
two-sided background estimation. GRB 090305A and 
120403A did not trigger on-board, and do not show 
compelling evidence for a signal in the offline data. 
GRB 090305A and 120403A also did not trigger on- 
board, but do show clear evidence of a signal present 
in the offline analysis (figure 5), though with relatively 
poor statistics. The remaining Swift short-GRBs from 
the observable sample both triggered on-board and are 
clearly identified offline. 
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Fermi sky at 090305 05:19:50.844 +0.512S 



spacecraft phi 


FIG. 5: Short GRB 090305A was seen by Swift and did 
not trigger GBM on-board. Low photon statistics for this 
weak event results in large sky location uncertainty shown 
by the 1, 2, and 3 a confidence regions. Arrows represent 
movement of the Fermi spacecraft over ±2 minutes. 


VI. CONCLUSION 

Direct detection of gravitational-waves from the 
merger of NS-NS or NS-BH binary systems is ex- 


pected to occur within the next few years as advanced 
2nd generation ground-based gravitational-wave de- 
tectors come online. Fermi GBM provides a unique 
and promising opportunity to observe EM counter- 
parts due to its large sky coverage, and the anticipated 
association between NS-NS mergers and sGRBs. To 
maximize the use of information from both GWs and 
GBM, we have outlined a joint search strategy trig- 
gered by observations of GW signals from coalescing 
binary systems in which a small amount of local GBM 
data is scanned using a likelihood-ratio based analy- 
sis applied to the full instrument data. We anticipate 
such an effort will allow sensitive follow-up of NS-NS 
and NS-BH mergers, most of which are expected to 
not be accompanied by triggered GRB observation. 
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